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LIFT DEVELOPED ON UNRESTRAINED RECTANGULAR WINGS ENTERING GUSTS AT SUBSONIC 

AND SUPERSONIC SPEEDS ^ 


By Harvard Lomax 


SUMMARY 

The object of this report is to provide an estimate, based on 
theoretical calculations, of the forces induced on a wing that is 
flying at a constant forward speed and suddenly enters a ver- 
tical gust. The calculations illustrate the effects of Mach 
number (from 0 to 2) and aspect ratio (2 to oo ) , and solutions 
are given by means of which the response to gusts having ar- 
bitrary distributions of velocity can be calculated. The ejects 
of pitching and wing bending are neglected and only wings of 
rectangular plan form are considered. Specific results are 
pi'esented for sharp-edged and triangular gusts and various 
wing-air density ratios. 

INTRODUCTION 

Studies of the gust-response problem for restrained wings 
(wings of infinite mass) entering sharp-edged gusts at super- 
sonic speeds are already well advanced. Miles, Strang, 
Biot, and Hcaslet and Lomax (refs. 1, 2, 3, and 4) presented 
solutions to such problems for two-dimensional wings; 
Miles and Goodman (refs. 5 and 6) presented solutions for 
rectangular wings having tip Mach cones that do not intersect 
the opposite edge. Miles and Strang (refs. 7 and 8) gave 
results for a triangular wing with supersonic edges. Theo- 
retical studies restricted to incompressible flow fields contain 
the classical solutions due to Wagner (ref. 9), Kussner (ref. 
10), and von Kdrmdn and Sears (ref. 11), the former con- 
taining the solution for the indicial lift on a two-dimensional 
sinking wing and the latter two containing the solution for 
the lift on a restrained two-dimensional wing entering a 
sharp-edged gust. The extension of these studies to include 
the gust response for wings of finite aspect ratio has been 
carried out by Jones (ref. 12). Later, further extensions to 
include the effects of gust shape as well as aspect ratio have 
been made by Zbrozek (ref. 13) and Bisplinghoff, Isakson, 
and O’Brien (ref. 14). 

The purpose of the present report is twofold: first, to 
present solutions for a two-dimensional restrained wing 
entering a sharp-edged gust at sonic and subsonic Mach 
numbers (specifically, Mach numbers equal to 1.0, 0.8, and 
0.5); and second, to use these results together with those 
mentioned above to estimate the effect of wing aspect ratio 
and airplane mass on the lift response for airplanes flying at 
various speeds through the Mach number range from 0 to 2 
and penetrating both triangular and sharp-edged gusts. 

A list of symbols is given in the appendix. 


ANALYSIS 

THE EQUATION OF MOTION 

If induced pitching moments are neglected, the motion of 
a rigid wing disturbed from its equilibrium position by 
arbitrary external lifting forces is governed by Newton’s 
second law. Thus, if w is the vertical velocity of the wing 
and m its mass, we can write 

m ^ = S forces (1) 

where the forces to be summed are the aerodynamic ones due 
to the gust velocity and the motion of the wing from its 
position of equilibrium. 

First, consider the force that results from a small vertical 
motion of the wing. Suppose the wing has been flying in 
steady level flight at a constant speed Ua up to a time f' = 0. 
Fix an xyz coordinate system in space (2 positive upward) 
such that at t' = 0 the y axis lies along, and the origin on, the 
wing leading edge and, further, such that the wing is moving 
in the negative x direction. For f'^0 the wing moves 
away from these coordinates, continuing forward at the 
constant speed along the negative x axis, and now also 
moving downward at a constant rate — w=Uoa. The 
transient lifting force on the wing induced by such a maneuver 
shall be referred to as the indicial lift (positive upward) and 
designated in coefficient form by the symbols Ci^ or foi' 
section or total lift values, respectively. 

Given the indicial lift coefficient, one can show by using 
the prinicples of superposition that the lift due to an arbitrary 
variation of angle of attack caused by the vertical velocity 
of the wing can be determined from the relation 

L(t')=-I^-^ £' C^^(t'-h')w(tf)dtf ' (2a) 

Next, consider the force that is induced on a wing pene- 
trating a sharp-edged gust having a uniform upward velocity 
Wo. If the wing is restrained so that it can move neither 
upward nor downward (corresponding in flight to the limiting 
case of infinite wing mass), the section or total lift coefficients 
induced by a unit value of WqIUq will be designated or 
Ol^, respectively. Figure 1 illustrates the differences be- 
tween the boundary conditions for, and the initial variations 
of Ci^ and Ci for a two-dimensional wing traveling at a 
subsonic speed. 


1 Supersedes NAG A TN 2925, “Lift Developed on Unrestrained Rectangular Wings Entering Gusts at Subsonic and Supersonic Speeds" by Harvard Lomax, 1953, 
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Figure 1. — Boundary oonditious and initial variations of ci^ and ci^. 

Given a value of Cz,^, one obtains the lift on a restrained 
wing flying into a gust having an arbitrary vertical velocity 
distribution, Wa{t'), by the relation 

Cz^it'-inwa(t,')dt/ (2b) 

Substituting equations (2a) and (2b) into equation (1), 
one finds the expression for the vertical motion of an unre- 
strained wing flying into a gust ; thus, 

If i ^?r,(i'-h')«^«(h')dh' (3a) 


Since WaiU') is assumed to be given, equation (3a) is an 
integral equation — in terms of w{t ') — of the second kind 
with a variable upper limit. It is convenient first to study 
equation (3a) when the gust is a step function (sharp-edged 
gust). For this case Wa{t\) becomes a constant Wq, say, and 
equation (3a) reduces to 



r' 

Uo dt' Jo 




qSwo 

Uo 


Cz,,in (3b) 


The solution to equation (3b) can be used to find the 
induced force on an unrestrained wing entering a gust of 
arbitrary structure. Methods for solving the integral 
equation and applying its solution will be developed in the 
subsequent sections. 

INDICIAL LIFT ON A SINKING WING 


The analysis involved in calculating the indicial lift force 
on the wing is based on the assumptions that underlie linear- 


ized, thin-au'foil wing theory in general. Mathematically, 
these assumptions imply that the governing partial differ- 
ential equation of the flow field is the three-dimensional wave 
equation. In terms of the velocity potential and for an 
axial system fixed relative to the still air at infinity, the wave 
equation can be written 


<Pxx-\-Vn~iripzi—<Pit (4) 

where t is the product of the speed of sound Uq and the time t'. 
For a wing moving in the s= 0 plane, the boundary conditions 
are that <p is continuous everywhere except across the wing 
and its vortex wake, is a constant over the region 

bounded by the wing plan form at any given time, and all 
velocities vanish outside the starting wave envelope. 

All values of and used herein have been presented 
in previously published reports. The indicial lift on a 
sinking rectangular wing traveling at supersonic speeds has 
been presented by Miles (ref. 5). As the aspect ratio tends 
to infinity, this solution approaches that for the two-dimen- 
sional case given in references 1, 2, 3, and 4. At sonic and 
subsonic speeds, results for a two-dimensional wing are 
available for Mach numbers equal to 1.0, 0.8, and 0.5 (ref. 15) 
and for the incompressible case (ref. 16). Finally, the 
indicial lift on sinking wings of finite aspect ratio in incom- 
pressible flow is presented in reference 12. Curves showing 
the effect of Mach number on the two-dimensional values 
are presented in figure 2 (a), and the effect of aspect ratio at 
supersonic speeds is indicated in figure 2 (b). Tabular 
values for the two-dimensional wing flying at Mach numbers 
equal to 0.5 and 0.8 are given in tables I and II, respectively. 



(a) Two-dimensional wing for several Mach numbers. 
Figure 2. — Variation of indicial lift response with chord lengths 

traveled. 
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0 12 3 4 

Chord lengths traveled, t'/c 

(b) Rectangular wings of aspect ratios 2, 5, and <». 

Figure 2. — Concluded. 


and the boundary conditions for a section in the xt plane are 
indicated in figure 3. 


Trace of foremost 



The solutions obtained for the load coefficient over regions 
1 and 2 shown in figure 3 can be written (details of the 
analysis are omitted): 

For region 1 


8 j Mo{t—x) /gj^) 

qg { Wq\ x(1 + Mo)V x + Mof 


For region 2 


RESPONSE OF A RESTRAINED WING TO A SHARP-EDGED GUST 

Load distribution, — The lift induced on a restrained wing 
penetrating a sharp-edged gust can also be determined by 
solving equation (4) subject to the proper boundary condi- 
tions.^ For a wing moving in the 3=0 plane, these condi- 
tions are similar to those given for the indicial lift on a sinking 
wing in that all velocities vanish outside the starting wave 
envelope and <p is continuous everywhere e.xcept across the 
wing and its vortex wake, but differ from the indicial case 
in that (<^ 2 ) 2=0 is a constant only over the portion of the 
wing plan form that bas penetrated the gust, being zero 
over the remaining portion (see fig. 1). This problem has 
been solved for a rectangular wing traveling at supersonic 
speeds by Miles (ref. 5) and, again, as the aspect ratio tends 
to infinity, this solution approaches that for the two-dimen- 
sional case given in references 1, 2, 3, and 4. Two-dimen- 
sional wings flying at ..he speed of sound and two- and three- 
dimensional wings flying in an incompressible medium have 
also been considered (refs. 4, 16, and 12, respectively). 

The problem of finding the two-dimensional gust response 
at subsonic speeds can be solved by the same method that 
was used in reference 15 to find the two-dimensional subsonic 
indicial response. For these cases equation (4) reduces to 

V’xx-hv>sz=Pii (5) 

3 It Is interesting to note that the gust lift function ci^ can be related to the indicial response 
following a step variation of angle of attack under quite general conditions by the reciprocal 
theorems given in reference 17. 


Ap^ 8 ( j Mo{t—x) 

qg w{l M o) \y x-\-Md 


2 

—^Mo(t — x){c — X — Mot) 

TV 


r 2K 

EF'(,p) +KE'(^k) -KF'{,p)-\ 

Lv(«^-x^(1-Mo^) 



(6b) 


The symbols E, K, E'{4'), and are elliptic integrals 

defined in the appendix, their modulus k being given by 


k 



(f a)(l -|- Mq) — 2c 
(f+:s)(l -{-Mo) 


and their argument by 

i/'=arc sin 




-|- Mot 


( 7 ) 

( 8 ) 


Equations (6a) and (6b) give the loading over the complete 
wing section for values of t less than or equal to 2cl{\—Mo). 
Since Motje equals Uot'lc, the number of chord lengths trav- 
eled, these equations represent the exact linearized solution 
for the section load distribution during the time required 
for the wing to travel 2Mo/(l-Mo^) chord lengths after reach- 
ing the gust front. Hence, for a Mach number equal to 0.8, 
equations (6) establish the gust response during the first 4.44 
chord lengths of travel after the gust penetration. Figure 4 
shows the variation of gust loading Apjqg throughout this 
interval and also, for comparative purposes, the indicial load 
variation Apiqa for the first four chord lengths of travel. 
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Figure 4. — Gust and indicial loading on two-dimensional wing. 

i\/o=0.8. 


The dashed curves in figure 4 represent the final steady-state 
load distribution adjusted so as to give the same total lift as 
the exact solutions for the gust and indicial cases at U4' Ic 
equal to 4.44 and 4, respectively. Thus, to the degi-ee of 
accuracy indicated in figure 4, the gust and indicial loadings 
at Mo=0.8 can be approximated for larger values of Ugt'lc 
by the expressions 

(9.) 


and 


> 

1 

4^ 

jc— X 

CiJJJ^'jc) 

qoi 

X 

- ^g(“) - 


d/o / Uoi' 
1-Mo^ c 


(9b) 


The variation of Ci^ and C;^ for values of Unt'lc greater than 
4 ivill be discussed presently. 

For a Mach number equal to 0.5, equations (6) are suffi- 
cient to establish the gust response for only the first 1.33 
chord lengths traveled. Further calculations were carried 
out and the exact loading was established for both the gust 
and indicial cases for values of Uat' jc less than or equal to 
2.33. These calculations were for the most part numerical 
and no simple closed expressions such as those presented in 
equations (6) were obtained. Figure 5 contains the I’esults. 
Again, the dashed curves represent the final adjusted steady- 


state load distribution indicating that the gust and indicia 
loadings for ilfo=0.5 can also be approximated for larger 
values of Unt'lc by the equations (9). 

Section lift. — Wlien integrated across the chord at a fixed 
time, the loadings shown in figures 4 and 5 give the variation 
of the lifting force on the wing section during the early portion 
of the response. In the interval Q<Uat' lc<Adoj{l-\-AdQ) 
equation (6a) integrates to give 



In the interval Ic<2MqI{\~A'Io) the c.x- 

pression for the loading is too complicated to integrate 
analytically and the section lift was calculated bj'^ numerical 
methods. The results, together with those for Ci^ (taken 
from ref. 15), are given in figure 6. Since, as time goes on, 
Ci^ must approach C;^, the emve for the gust response was 
simply faired into the curve for the indicial response in the 
manner shown by the dashed lines. Finall 3 '‘^ for values of 
Uot'lc greater than 10, the following equations, taken from 
reference 1, can be used; 

For Ada— 0.5, 


27t (, 1.333 44.218 ) 

( 5-f207of7c) {5 + 2{Uot'lc)f] 


(11a) 




0 2 4 6 8 10 

Chord lengths traveled, UJ'/c 

(a) Two-dimensional wing for several Mach numbers. 
Figurb 7. — Response of restrained wing to unit sharp-edged gust. 
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(b) Rectangular wings of aspect ratios 2, 5, and <» . 
Figure 7. — Concluded. 


The final curves (determined from the previous analysis 
and aforementioned references) for Ci^, the section lift coeffi- 
cient developed on a restrained wing entering a sharp-edged 
gust, are shown in figure 7(a) for Mach numbers equal to 
0, 0.5, 0.8, 1.0, 1.2, 1.41, and 2.0. Tabular values are given 


in tables I and II for Mach numbers equal to 0.5 and 0.8. 
The effect of aspect ratio at a Mach number equal to 1.41 
is shown in figure 7(b). 

TABLE II.— VALUES OF ci„ AND c, FOR Mo=0.8 


Udt'Ic 




Vl-Afo2ci,/2. 

0 

0. 478 

0 

.1 

.466 

.044 

.2 

.454 

.085 

.3 

.442 

.129 

.4 

.430 

.170 

.5 

.423 

.209 

.6 

.426 

.234 

.7 

.433 

.256 

.8 

.442 

.276 

.9 

.451 

.296 

1.0 

.461 

.315 

1.5 

.507 

.402 

2.0 

.546 

.465 

2.5 

.581 

.513 

3.0 

.610 

.551 

3.5 

.632 

.584 

4.0 

.662 

.616 

4.5 

.670 

.642 

5.0 

.687 

.663 

6.0 

.714 

.700 

7.0 

.738 

.730 

8.0 

.760 

.758 

9.0 

.779 

.780 

10.0 

.798 

.796 

CO 

1.000 

1. 000 


RESPONSE OF AN UNRESTRAINED WING TO A SHARP-EDGED GUST 

Given the indicial lift response and the response for a 
restrained wing penetrating a sharp-edged gust Ci,^, one can 
use equation (3) to find the motion of an unrestrained wing 
entering a sharp-edged gust having a constant upward 
velocity Wq. As in reference 4, the lift on the unrestrained 
wing can be related to an infinite series of integrals involving 
and First, set 

Uot' _ Uoh' 

w 2m 

“77o““’ ^^pocS 

SO that, by integrating equation (3) with respect to t' , one 
finds 


TABLE L— VALUES OF c,„ AND c,^ FOR Mo=0.5 


Uot'/c 



Vl-Afo!o„/2, 


0 

1. 103 

0 

.1 

.995 

.081 

jy 

.882 

.158 

.3 

.772 

.236 

.4 

.691 

.301 

.5 

.656 

.347 

.6 

.634 

.382 

. 7 

.618 

.413 

.8 

.604 

.441 

.9 

.588 

.438 

1.0 

.579 

.488 

1.1 

.581 

.509 

1.2 

.597 

.526 

1.3 

.617 

.545 

1.4 

.634 

.560 

1.6 

.642 

.673 

1.6 

.659 

.586 

1.7 

.669 

.598 

1.8 

.677 

.609 

1.9 

.685 

.619 

2.0 

.693 

.629 

2.5 

.724 

.673 

3.0 

.746 

.704 

3.5 

.765 

.727 

4.0 

.781 

.746 

4.5 

.795 

.764 

5.0 

.806 

.779 

5.5 

.816 

.793 

6.0 

.825 

.806 

6.5 

.832 

.815 

7.0 

.840 

.826 

7.5 

.848 

.835 

8.0 

.854 

.845 

9.0 

.866 

.862 

10.0 

.877 

.877 

CD 

1.000 

1.000 


■U!=— r Gi (Ti)dri-| r Ci, {t Tl)W{Tl)dTi 0 (12) 

M Jo * Mdo 

Then use the relation 

■ ^,if = {dw/d t) in I Wo) 

WolUo 

and iterate equation (12) using Liouville’s method of suc- 
cessive substitutions. (See ref. 18.) This yields 


-%7==C'i (r)— - r — Ti)(7i (Ti)dn-4- 

Wo/Uo * mJo “ * 

\ (t — Ti)dri r ' — T2)(7i^(r2)dT2— . . . (13) 

M Jo “ do 

Equation (13) converges uniformly ^ for all t. By means 


3 The statement made in reference 4 on the convergence of this series is unnecessarily 
restrictive. Since the greatest values of (7i,„ and are (=>) and Cl, (»), the general term 

of the series does not exceed 




” ■ • -Jo 


that is, does not exceed 

i[Cr.„(«)]"Ci,(») At- 
and by the ratio test the series converges uniformly. 
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0 2 4 ,6 8 10 12 14 

Chord lengths traveled, Ugt'/c 


(a) Mo=0, A = “ 

Fiourb 8. — Response of unrestrained wing to a uniform, sharp-edged 

gust. 



(d) Mo=\.G. A=a> 

Figure S. — Continued. 



(b) Mo = 0.5, yl = 
Figure 8. — Continued. 



0 2 4 6 8 10 12 14 

Chord lengths traveled, U^t'/c 



(e) Mo=1.2, A=oo 

Figure 8. — Continued. 



293052—51 2 


(c) Mq=0.8, A=<» 
Figure 8. — Continued. 


(f) il/o= 1.41, .4=0° 
Figure 8. — Continued. 
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(g) A/o=1.00, ^ = 5 
Fkujre 8. — Continued. 



(i) M„=1.41, A = 5 
Figure 8. — Continued. 



(h) il'/o=1.02, yl = 5 
Figure 8. — Continued. 



(i) M„ = 1.2, A = 5 
Figure 8. — Continued. 



(k) d/o = 1.12, A = 2 
Figure S. — Continued. 



(1) 71fo=1.41, A = 2 
Figure 8. — Concluded. 
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of it, CJ{wqIUo) and Ci/{wolUo), the total and section lift 
coefficients induced on unrestrained wings entering a sharp- 
edged gust, have been calculated and the results are shown in 
figure 8. Table III indicates the range of Mach numbers and 
aspect ratios for which calculations were made, the numbers 
in the chart referring to the individual figures in which the 
results are presented. It should be noted that results are 
given for a wing flying at Ma=l and having a finite aspect 
ratio. Such cases can be calculated from the indicial and 
restrained gust responses presented in reference 5. These 
responses are still valid at iIfo=l for values of the time 
variable up to that for which the wave envelope induced by 
one side of the wing crosses the opposite side. A wing of 
aspect ratio 5 flying at the speed of sound travels 13 chord 
lengths during this time interval, and this is sufficient to 
establish the significant part of the response to a sharp- 
edged gust for /ii<300. 

TABLE III.— VALUES OF MACH NUMBER AND ASPECT 
RATIO FOR WHICH CALCULATIONS WERE CARRIED OUT 


AA \ 

\ 

2 

5 

CO 

0 



8a 

. 5 



8b 

.8 




1.0 


8g 

8h 

8d 

1.02 


1. 12 

8k 


1. 20 

8i 

8e 

1.41 

81 

8j 

8f 


The chart also shows that the gust response for the 
unrestrained wing was calculated at Ma~0 for an infinite- 
aspeot-ratio wing (for comparative purposes) but not for 
finite-aspect-ratio wings. The gust response onjboth infinite- 
andi finite-aspect-ratio wings in incompressible flow has been 
studied extensively by means of operational methods in 
references 12, 13, and 14. Wliere comparisons can be made, 
the results obtained in this report using'equation (13) agree 
well with those given in the references mentioned. 

RESPONSE OF UNRESTRAINED WING TO ARBITRARY GUST 

The function ClKwoIUo) presented in the previous section 
can be thought of as the indicial gust response for lift on an 
unrestrained wing. In this sense it is apparent that the 
lift on a wing penetrating a gust in which w is a function of 
the chord lengths traveled can be calculated by superposi- 
tion and is represented by the integral 


Gl — 


A 

dr 


Jo WqIUo Uo 


dn 


(14) 


By means of equation (14), the lift induced on a wing moving 
at the constant speed Uo and entering a gust, the vertical 
velocity of wliich starts at zero and increases linearly with 
distance of penetration, is simply the integral of CJ{woIUq). 
Thus, representing the section lift coefficient developed by a 
wing entering a gust with a unit gradient by the symbol Ci,^, 
we can write 



dri 


( 15 ) 


If the vung flies into a gust with a triangular-shaped dis- 
tribution of w, having its maximum intensitj’' Wt a distance h 


chord lengths from the front, it follows at once that the 
resulting lift response Cz.l(w,IUo) is given by 


Cl 

WtlUo' 


^\[CLS.r)]-0<r<h 

\ [OLj,r)~2C^JiT-h)]- h<r<2h 

\[CLS.r)-2CLS.r-h)^ 

\^CL<,r-2h)\,2h<T 


^ (16; 


Examples of the various gust shapes and the responses in 
lift and vertical motion of wings penetrating them are shown 
in figure 9. 


Variation of gust intensity ■ 
Path of point on wing 
Variation of lift on wing 




Figure 9. — Gust shapes and responses. 


MAXIMUM LIFT DUE TO GUST PENETRATION 
SHARP-EDGED GUST 

Consider the maximum increase in lift caused bj^ the 
entry of the wing into a sharp-edged gust. This increment 
is given for the range of Mach numbers, aspect ratios, and 
wing-air density ratios shown in table III bjr the maximum 
values of C^l (wolUo) and CiKwqIUo) on the curves shown in 
figure 8. 

First, let us consider wings of infinite aspect ratio. For 
such wings the variation of the maximum gust-induced lift 
coefficient with Mach number is shown in figure 10, and a 
cross plot in which Mo instead of ja is held constant is pre- 
sented in figure 11. The values for ;i= oo are the steady-state 
values given bj^ the simple equations 



2irl 1 — ^Mo'^ 1 
il-\JMo ^ — 1 ; M(G> 1 


The difference between the lift increment on a restrained 
wing and that on one with a finite value of ja is seen to be 
most pronoimced at the high subsonic M ach numbers. Notice, 
for example, that the percentage increase in [Ci/{'Wo/Uo)]max 
found by increasiag from 0 to 0.8 is 67 for «=oo 
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0 .2 4 .6 .8 , 1.0 1.2 1.4 1.6 

Mach number, 

Ficure 10. — Maximum increment of lift induced on a two-dimensional wing entering a uniform sharp-edged gust. 



0 100 200 300 


Wing -air density ratio, 

Figure 11. — Maximum gust-induced lift on two-dimensional wings. 

(Prandtl Glauert rule) but only 37 for p=200. Table IV 
indicates the relative increase in [ci/(wo/f^o)]mM caused by 
compressibility for three different values of the wing-air 
density ratio. 


TABLE IV.— PERCENT INCREASE IN (c,)„„3 RELATIVE TO 
ITS VALUE AT Mo = 0. 


\ 

X P- 
\ 

Mo \ 

\ 


300 

200 

0 

0 

0 

0 

.8 

67 

39 

37 

1.0 

e» 

103 

88 

1.2 

-4 

10 

12 


Consider next the effect of aspect ratio on the maximum 
lift increment induced on a rectangular wing penetrating a 
sharp-edged gust. When p= «> this increment is again given 
by the steady-state value of the lift-curve slope and is 
presented for A=°^, 5, and 2, in figure 12. These steady- 
state values are taken from the numerous studies made of 
lifting surfaces traveling at subsonic and supersonic speeds. 
On the supersonic side, for cases in which 1 >1, the 

equation (see, e. g., ref. 19) 

4 ^ 

\Wo/Uo/ max VMo"- 1 \ 2 A VMo'- 1 / 
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Figure 12. — Maximum increment of lift induced on restrained rectangular wing entering a uniform sharp-edged gust. 


applies and, for cases in which 1 >A-\lM^— \ >0, the curves 
in reference 20 were used. On the subsonic side the portions 
of the curves in the range Q<AyJl—Mo < 2 were again taken 
from reference 20. The results in references 19 and 20 are 
sufficient to cover the entire Mach number range for tlie 
A = 2 wing. For the A = 5 wing the values on the subsonic 
side outside the range 0 <2 were taken from a 

curve * that was compiled from a large number of solutions 
for lifting surfaces traveling at subsonic speeds. 

The values of for rectangular wings 

traveling at supersonic speeds, given in figure 8, and the 
incompressible-flow solutions, given in references 12, 13, 
and 14, were used to prepare the curves in figure 13. The 
dashed lines between the Mach numbers of 0 and that for 
which A-^Ma—l = l are interpolated, the two-dimensional 
results presented in figure 10 lending credence to the validity 
of the interpolation. Figure 14 presents the aspect-ratio 
effect on [CJiwolUAlmax at Mo==l. 

It should be noted that in the vicinity of Mo= 1 , the curves 
for which A= o6, fj,= oo (figs. 10 and 12), and probably also 

< The curve was taken from an article prepared by Robert T. Jones and Doris Cohen for 
the forthcoming series on High-Speed wALerodynamics and Jet Propulsion, Princeton Uni- 
versity Press. 


those for which A— 5, ti= <» (figs. 12 and 13 (a)), are not 
valid representations of the gust-induced lift on actual wings 
flying at these speeds, although they do represent solutions 
to equation (4) consistent with the boundary conditions 
previously discussed. For the two cases mentioned, the 
assumptions on which equation (4) is based are violated. 
These assumptions are more closely approached, however, 
as the wing-air density ratio and the aspect ratio decrease. 
Hence, for lower values of A and ^ the solutions given herein 
for wings traveling in the transonic speed range have justifi- 
cation on a physical as well as a mathematical basis. 

TKIANGULAB GUST 

The maximum inci’ease in lift on a two-dimensional wing 
passing through a triangular gust having its maximum 
intensity 12 chord lengths from its front is shown in figure 15. 
For the lower values of ix, the variation of [cr/(w,/Uo)]mai 
with Mach number is similar to that calculated for the 
sharp-edged gust and shown in figure 10. As ix increases, 
however, a comparison of the results shown in these two 
figures indicates the importance of the assumed gust shape in 
estimating the maximum gust-induced lift. Table V shows 
the difference in the compressibility effect obtained for the 
sharp-edged and triangular (12^chord lengths to apex) gusts. 
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(a) A = 5 

Figure 13. — Maximum increment of lift induced on a rectangular iviiig entering a uniform sharp-edged gust. 


TABLE V.— RATIO OF THE VALUE OF AT Mo=0.8 AXD 

1.0 TO ITS VALUE AT Mo=0. 


Gust 
\ type 

. \ 

\ 

\ 

Mo= 

=0.8 


= 1.0 

Sharp 

edge 

Tri- 

angle 

Sharp 

edge 

Tri- 

angle 

60 

1.24 

1.24 

1.44 

1.51 

100 

1.20 

1.29 

1.59 

1.58 

200 

1.37 

1.33 

1.88 

1.67 

3C0 

1.39 

1.34 

2.03 

1.70 

03 

1.66 

1.38 

03 

1.85 


Figure 16 presents the aspect-ratio effect on the maximum 
lift response for rectangular wings penetrating triangular- 
shaped gusts. The values at Mo=0 were calculated from 
the results given in reference 12 and again the dashed lines 
represent an interpolation. 

CONCLUDING REMARKS 

Results are presented for the lift developed by a restrained 
two-dimensional wing flying at a Mach number equal to 0.5 


or 0.8 and penetrating a sharp-edged gust. Similar results 
are reviewed for Mach numbers equal to 0, 1.0, 1.2, 1.41, 
and 2. 

A method is given whereby the lift can be estimated 
(neglecting the effects of airplane pitching and wing bending) 
for unrestrained rectangular wings in the aspect-ratio range 
2 to 00 , flying in the Mach number range 0 to 2, and penetrat- 
ing gusts of arbitrary structure. Specific results are given 
for sharp-edged and triangular-shaped gusts. 

In general, given variations in the wing aspect ratio, the 
wing-air density ratio, and the gust shape have their maxi- 
mum effect on the gust lift when the wing is flying at a 
high subsonic speed. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feh. S, 1953. 
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(b) A=2 

Figure 13. — Concluded. 


/^o 


max 



100 200 300 


Wing-air density ratio, 


Figure 14. — Maximum gust-induced lift on rectangular wings. 
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Fiquee 16. — Maximum increment of lift induced on a two-dimensional -nung entering a triangular gust having its apex 12 chord lengths from the 

front. 




Figure 16. — Maximum increi 


r/iMM —M— ll 
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APPENDIX 


LIST OF SYMBOLS 


aspect ratio 
speed of sound 
chord length 

wing lift coefficient, 


qS 


section lift coefficient, — 
qc 

complete elliptic integral of second kind with 
modulus k 

incomplete elliptic integral of second kind with 


modulus k' and argument i}/ 
incomplete elliptic integral of first kind with modu- 
lus k' and argument \j/ 


Wq 

Uo 

number of chord lengths from front to apex of 
triangular gust 

complete elliptic integral of first kind with modu- 
lus k 

modulus of elliptic integrals (See eq. (7).) 
lift on wing 

Mach number at which wing is traveling 
mass of wing 

loading coefficient, pressure on lower wing surface 
minus pressure on upper wing surface divided 
by dynamic pressure 

djmamic pressure, ^ poW 

wing area 

time 

aot' 

horizontal velocity of wing 
vertical velocity of wing 
velocitj^ of arbitrary gust 
velocit}' of uniform, sharp-edged gust 
maximum velocitj’ of triangular gust 
Cartesian coordinates fixed with reference to still 
air at infinity, 2 positive upward, y parallel to wing 
leading edge, negative x direction corresponding 
to direction of wing motion 
wing angle of attack 

2 jjfi 

wing-air density ratio, ^ in analj^sis of complete 

PoCA 

2 Tfh 

wing, — 2 in analj^sis of wing section 

PoC 

air density 

U t' 

chord lengths traveled by wing, — ^ 

c 

perturbation velocity potential 
argument of elliptic integrals (See eq. (8).) 


Subscripts: 

g response of restrained wing to unit, sharp-edged 

gust 

s response to gust with unit velocity gradient 

a indicial response on sinking wing 
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